Incentives and regulations in bike-sharing 
systems with stations of finite capacity 

Christine Fricker* Nicolas Gast^ 
January 6, 2012 



Abstract 

As of 2011, more than 200 cities around the world have a bike sharing 
system. In such a system, users arrive at a station, use the bike for a while 
and return the bike to another station. This paper presents a model taking 
account of the finite number of bike locations at the stations. In case of 
symmetry, as the system gets large, the mean field limit provides an insight 
of the system behavior. Convergence of the invariant measures is proved 
and closed form results are obtained. The influence of the parameters and 
of various load balancing strategies on the performance, measured by the 
proportion of so called problematic stations, is discussed. Even in this 
symmetry case, the system exhibits a poor performance. We show that 
simple incentives, such as suggesting users to return to the least loaded 
station among two, improve dramatically the situation. An asymmetric 
scenario is also investigated. In that case, simple incentives are not enough 
and regulation mechanisms, such as redistribution of the bikes by trucks, 
are needed. 

1 Introduction and related work 

Vehicle rental networks are giving rise to interest in the last few years. Shared 
car networks were hard to develop, but bike sharing networks have been deployed 
successfully in many cities giving rise to different research activities. 

History 

The first bike-sharing system started in July 28, in 1965 in Amsterdam with the 
Witte Fietsen but the bikes where not locked and the system collapsed within 
days [2] . Few other systems have been developed until Paris launched the Velib 
program in July 2007, which helped popularize the idea. Velib consists of 20 
000 available bikes that are hired at 1 500 stations. The stations can be freely 
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accessed at any time and the bikes returned back at any station. There are now 
such other programs in around 200 cities like Velo'v in Lyon since May 2005, 
4 000 bikes at 334 stations, Bicing in Barcelona since March 2008, Vel'oh in 
Luxembourg, Bixi in Montreal with now 5K bikes at 400 stations or Ecobi in 
Mexico. Now the greatest one is in Hangzhou, China, started in May 2008, with 
50K bikes at 2000 stations, one every 100 meters. 

Description and related work 

In a vehicle rental network, users arrive at rental stations, use a vehicle and then 
return the vehicle to any rental station. Several papers have studied vehicles- 
sharing system focussing on economical aspects, growth trend or environmental 
impact, see for example [7, 11, 3, 2] and the references therein. However, few 
analytical models are proposed so far. [1] developed a simulation model to study 
the impact of various parameters such as the fleet size or reallocation policies. A 
closed queueing network is also proposed which has a product form steady-state 
distribution but studying the system getting large is difficult because of the 
combinatorial explosion of the expression of the normalization constant which 
is a partition function. An approach is to compute an asymptotic expansion of 
this partition function when both the number of stations TV and customers M 
get large with the ratio M/N tends to some constant s (cf [G] for a review of 
the literature). 

A different point of view is adopted George and Xia in [ ], remaining the 
number of stations fixed and letting the fleet size get large. They obtain, with 
much elementary technique, the behavior in the subcritical regime in some sense, 
j )] by Malyshev and Yakovlev deals with large Jackson networks. Their result 
is that there exists a critical value s cr for the fleet size to number of stations 
ratio s. If s < s cr then the mean number of customers in the stations is finite. If 
s > s cr then some stations are bottlenecks concentrating a number of customers 
which tends to infinity while, in the other stations, the number of customers 
remains with finite mean and a geometric distribution. These results are ob- 
tained analytically using the product form expression of the invariant measure 
and give a good idea of the behavior of the system, see also [ ] for a probabilistic 
approach. 

Nevertheless, these models do not catch one of the main features of the 
application, the fact that the stations have a finite number of bike locations. 
For the user, it means that, when he returns the bike, the station could be 
saturated. The main purpose of the paper is to build an analytical model and 
investigate the case where the stations at which the bike can be hosted have 
finite capacities. 

Contribution 

We consider a model of a bike-sharing system and we study the steady state 
performance of this system. The system is composed by a large number N 
of stations and a fleet of s ■ N bikes. Each station can host up to a finite 
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number of bikes K called its capacity with network terminology throughout of 
the paper. We study both a symmetric case where all stations are identical and 
an asymmetric case where some stations are more popular than others. We study 
various incentives or regulation mechanisms that can improve the performance. 
Since no closed form expression for the invariant measure is available for the 
system with a fixed size, we use mean field techniques to prove that the limiting 
invariant measure as the system gets large can be computed easily as the Dirac 
mass at the equilibrium point of a differential equation. 

We first study a simple model with symmetry where every stations have 
the same parameters. We prove the existence of a Lyapunov function for the 
limiting dynamical system which helps to obtain the convergence of the invariant 
distribution to a Dirac mass. A critical value for the limit number of bikes 
per stations s is obtained, which has a simple interpretation for the limiting 
behavior of the system, and which corresponds to the optimal fleet size for our 
performance metric. This value is given in a closed form as a function of the 
capacity K of the stations and the arrival rate A at each station and the average 
trip duration n~ l . The optimal fleet size is K/2 + A//i and the proportion of 
problematic stations is 2/(K + 1) which decreases slowly with K. 

To improve the situation, we study different policies in case of saturation of 
the stations. Two directions are investigated: Firstly, in practice, technology 
allows users to access to real time information as the empty and saturated sta- 
tions, or more generally the state of the stations. The user could have incentives 
to choose the station where he returns the bike, maybe due to pricing. Secondly 
regulation is introduced, by trucks added to move the bikes to balance the sta- 
tions. It is mainly prove that returning the bike to a non saturated station 
does not change significantly the behavior of the stations and the performance 
with our metric, but returning the bike at the least loaded station among two 
improves the performance, even if just a fraction of the users do this. For the 
regulation, it depends a lot on the number of trucks per station that are de- 
ployed. Mean field limits are again used. Results give the limiting invariant 
measure in the first case as the system gets large. 

In order to break symmetry, a model with clusters is proposed, each cluster 
representing a level of popularity of a class of stations. We study the perfor- 
mance of regulation techniques and incentives by a numerical evaluation of the 
limiting dynamical system. These results show that having a mechanism to bal- 
ance the load between popular and non-popular stations is necessary to obtain 
a good performance. 

Organization of the paper 

Section 2 deals with the model description and the mean field results. In Sec- 
tion 3, the basic model with symmetric stations with finite capacity is completely 
analyzed. We study incentive mechanism in Section 4, where a fraction of users 
choose the least loaded of two stations to return their bike. Section 5 investi- 
gates a regulation where bikes are moved by trucks. This regulation is compared 
with the incentive mechanism, both in a symmetric system and an asymmetric 
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scenario. To improve the readability of the paper, most of the technical proofs 
are delegated to the appendix A. 

2 Mean field limits, stability 

This section presents the basic simple model which is studied. It is described 
with clusters even if most of the paper deals with one cluster. It helps to 
understand the target model for the applications. This model admits a mean 
field limit. It is formalized via the convergence of the empirical measure process 
to some dynamical system. It gives a strong framework to study the large scale 
bike systems: the concentration of invariant measures is shown, avoiding the 
presence of oscillations of the system. Thus studying the limiting behavior in 
Section 3 makes sense. The steps to obtain the concentration result are detailed 
here. All the other examples will fit roughly the same frame, but sometimes 
achieving the result fails because of technical gaps. 

2.1 Bike sharing model 

We consider a Markovian model of a bike-sharing system with N stations and a 
fleet of s ■ N bikes (s bikes per station in average). A bike can be either hosted 
at a station or in transit between two stations. Each station can host up to K 
bikes. The stations are grouped in C > 1 clusters, each cluster can correspond 
to a location in the city or to a level of popularity of some stations. There are 
N c stations in cluster c and all stations in the same clusters share the same 
characteristics. Assume when N gets large that N c /N tends to some constant 
a c . 

At each station of cluster c, new users arrive at rate A c . If there is no bike 
at the station, the new user leaves the system unhappy and chooses another 
mode of transportation. If the station is non empty, the user takes a bike at 
this station, chooses another station at random and rides to this station. The 
probability for this station to be in cluster d is p cc i and the station is choosen 
at random inside a cluster. 

The travel time from a station of cluster c to a station of cluster d is expo- 
nentially distributed of parameter \i cc i . A user arriving at destination wants to 
return his bike. If the station has less than K bikes, the user returns his bike 
and leaves the system. If the station has K bikes, no more bike can be hosted at 
this station (we say that this station is saturated). In that case, the user keeps 
his bike and chooses another station in cluster d at random and rides to this 
station with a travel time with an exponential distribution with parameter fj, cc . 
This process is repeated until he finds a non-saturated station. 

2.2 Mean field limit and steady-state behavior 

In this section, we show that the behavior of the system greatly simplifies as N 
grows. In order to simplify the notations, we restrict our presentation to the 
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one cluster case. It can be adapted straightforwardlly in the multiple cluster 
case, at the price of more complicated equations. 

As we consider a system with only one cluster c, for simplicity, we denote 
A = A c and /i = fi cc . Due to symmetry with respect to the N stations, let us 
define the proportion of stations with k bikes by 



Y k if) = (number of stations with k bikes at time t) (1) 




Because of the symmetry assumption (stations are chosen at random within a 
cluster), the process (Y (t)) — (Y^ (i),0 < k < K) is a Markov jump process 
with state space 

y N = {y = iy^--.,y K )e(^ 

The set y N is included in the set V({0, . . . ,K}) of probability measures on 
{0, . . . , K}. The condition J2ke{o k} kykN < sN means that the number of 
bikes at the stations is less than the total number of bikes in the system is 
sN. Therefore, there are sN — J2ke{o k} kykN > bikes that are in transit 
between two stations. 

If for k E {0, . . . , K} there is a proportion y k of stations with k bikes in the 
system, the arrival rate of users in a station with k bikes is XNy k . If k > 1, 
such an arrival causes yk to be decreased by 1/N and yk-i to be increased by 
1/iV. Similarly, if for each k there is a proportion yk of stations with k bikes, 
there are N{s — X)fce{o k} ky k ) bikes in transit between two stations. Since 
every bike arrives at rate /z, there is an arrival in a station with k bikes at rate 
VkfJ-Nis - T,ke{a,...,K} k Vk)- 

Denote by the fc-th unit vector of the transitions of Y N (t) can be 

summarized by: 

y -> y+ ~^i e k-i - e k ) at rate Xy k Nl k>0 



N 

y ->• y + jTfi e k+i - efc) at rate yk^jsN - ^ ny n N)l k<K . 



K 



Although the behavior of this system is complicated, Y N is a density de- 
pendent population process (see Kurtz [8]): there is a (finite) set C C 1SL K such 
that each transition y —> y + l/N occurs at rate N^iiy). This implies that 
the average change in a small interval dt is f(y)dt — J2eec £Pi{y)dt- Using the 
framework of Kurtz [ ] and the fact that / is Lipschitz-continuous, it proves that 
as N goes to infinity, for each T > 0, the process (Y N (t), < t < T) converges 
in distribution to (y(i),0 < t < T) unique solution of the following differential 
system of equations: 

K 

»(*) = Yl ( A ( e ^i " e k)lk>o + H(* - i J > y(t)])( e k+i - e k )l k<K ) (2) 

k=0 
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where [I,y] = T,ke{o,...,K} k Vk- 

The first term corresponds to the rate of arrival of new users to a station 
and the second term to the rate users return the bikes, which is fi times the 
proportion of bikes in transit at time t. 

The above differential equation (2) on y rewrites 

m = y{t)L v( t) (3) 

where the second member of equation (3) is the product of the probability mea- 
sure y, a row vector, by the jump matrix L y . For y € y, L y is the infinitesimal 
generator of an M/M/l/K queue with arrival rate /i(s — [I, y]) and service rate 
A. 

Equation (3) contains the mean field property of the model: it means that, 
when N tends to infinity, the empirical distribution y(t) of the stations evolves in 
time as the distribution of some non-homogeneous Markov process on {0, . . . , K}, 
whose jumps are given by L y u\, updated by the current distribution y(t). These 
jumps rates are those of a M/M/l/K queue, where the arrival rate fi(s— [I, y{t)]) 
is time dependent and the service rate is A. This queue represents the instanta- 
neous evolution of a "typical" node, because all the nodes have the same evolu- 
tion due to the symmetry. So for N large, this "typical" node reports the whole 
network. All the one cluster models that follow will fall into this framework, 
except the model with the choice which will be constructed slightly differently. 

Throughout the paper, we investigate the steady-state behavior of the sys- 
tem. A common way to do this is to compute the equilibrium point of the 
ODE (2), which is often unique, and then, in this case, as N grows, to prove 
that the steady-state distribution of process (Y (t))t>o concentrates around 
this unique equilibrium point y. However, this result has to be handle with care 
since the ODE may have limit cycles or chaotic behavior. 

In fact, it is quite standard to prove that, if all trajectories of the ODE 
converges to a unique equilibrium point of the dynamical system, then the 
sequence of invariant measure of Y N - which exists as Y N is an irreducible 
Markov jump process on a finite state space y N - concentrates on this fixed 
point as N grows. It means commutation of the limits as N and t tend to 
infinity. In practice, it avoids long term oscillations of the system, which would 
prevent to guarantee a given performance. For that purpose, finding a Lyapunov 
function is a possible tool and not so easy. For the model, the generic form of 
the limiting rates will be crucial to exhibit a Lyapunov function. It will be 
developed in Section 3. 

2.3 Performance metric 

When a station is empty, a new user arriving at this station cannot take a bike 
and leaves the system. Similarly, when a station is saturated, a user that wants 
to return his bike has to find another station. In practical situation, a high 
proportion of saturated stations causes more problems than empty stations: 
While a user finding no bike can leave the system, a user that cannot find a 
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spot for his bike has to find a spot somewhere else. If we want people to rely 
on this transportation mechanism, both proportions should be very low. 

Therefore, in all of the paper, we mainly focus on one performance indicator, 
called the limiting proportion of problematic stations. For N fixed, the propor- 
tion of saturated stations plus the proportion of empty stations in steady state: 
Y Q N + Y$ where y^(oo) is denoted by Y N . We prove that the ODE (2) has a 
unique equilibrium point y. to which all trajectories converge. Thus the Dirac 
mass at y is the concentration point of the sequence of the invariant measures. 
Then, our performance metric can be expressed as 

lim Y N + Y$ = y + y k , 

N— foo 

where lim means the limit in distribution, which is the limiting proportion of 
problematic stations as the system gets large. 

3 The symmetric model analysis 

This section is devoted to the basic one cluster model. As seen in Section 2.2, 
the behavior of the system can be approximated by the ODE (2) when the 
system gets large. The interpretation of the ODE as equation (3) which reflects 
the evolution of a typical node is used to first prove that the ODE has a unique 
equilibrium point to which all trajectories converge. This shows that the Dirac 
mass at this point is the concentration point of the invariant measures of the 
system as the number of stations gets large. This allows a complete study of 
the performance metric as a function of s. An optimal fleet ratio is found and 
the best performance is derived. The influence of parameters X/fi and K will 
be discussed. 

3.1 Steady-state of the one cluster model 

By equation (3), an equilibrium point of the ODE (2) is a solution of yL y = 0, 
where L y is the infinitesimal generator of an M/M/l/K queue with arrival rate 
(jl(s — [I,y]) and service rate A. For p > 0, let v p be the invariant probability 
measure of a M/M/l/K queue with arrival-to-service rate ratio p. For p = 1, 
v p is the uniform distribution on {0, . . . , K} while, for p ^ 1, v p is geometric: 
v p (n) = P n {Z{p))- 1 for all < n < K where Z(p) = E„ e{ o,...,K} P" = i 1 ~ 
p K+1 )/(l — p) is the normalizing constant. 

Therefore, the equilibrium points of (3) are given as the solutions of the 
fixed point equation 

V = v p{vY ( 4 ) 

where p(y) = fi(s - E fce {o,...,if} k Vk)/^ 

The key point is to prove the existence and uniqueness of the equilibrium 
point. Equation (4) helps a lot. Since for each p, there exists a unique stationary 
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measure v p , this equation is equivalent to: 



^{s-[I,v p ])=p. (5) 

But the mean [/, v p \ of the distribution v p is strictly increasing in p implying that 
p i y p(s — [/, v p ])/\ is strictly decreasing on R+, positive at 0. Thus, for each 
s > 0, there is a unique p € solution of equation (5). Thus there is a unique 
equilibrium point v p for the dynamical system (3) where p is characterized by 
fixed point equation (5) which can be rewritten: 

s = -p+ [I,v p ]. 
M 

3.1.1 Convergence of the trajectories to the fixed point 

For y and y' two probability measures on {0, . . . , K}, let h(y/y') be the relative 
entropy of y with respect to y' defined as 

h(y/y')= £ y„log^. (6) 

ne{0,...,K} Vn 

Recall that L y is the infinitesimal generator of an M/M/l/K queue with arrival- 
to-service rate ratio p(y) = ip([I,y}) where ip is positive C 1 on [0, s] defined by 
ip(x) = p(s — x)/\. This expression shows that p{y) is a function of [I,y]- It 
leads to adapt Proposition 4.1 in [12] which gives the existence of a Lyapunov 
function for the dynamical system. The reversibility of the generator L y is 
crucial. Note that [12] deals with communication networks. 

Theorem 1. Let ip be the function x i— > J x<p' (x) / ip(x)dx and g be defined as 

g{y) = h{y/v p{y) ) - log Z{p{y)) + </>([/, y]). 
Then g is a Lyapunov function for the dynamical system: 

V* : -Jj.9(y(t)) < and ^9(y(t)) = implies y(t) = v p(y) . 

Proof. The proof is given in Appendix A.l. □ 

A consequence of this theorem together with the uniqueness of the fixed point 
V = v p(y) 1S > as recalled in Section 2, that the steady state distribution of the 
system concentrates on the Dirac measure on y, as the number of stations grows. 
In practice, this facilitates greatly the study of the steady-state performance. 

3.2 Proportion of problematic stations 

Let y be the equilibrium point of the equation. When the number of stations 
is large, the proportion of empty stations in the system can be approximated 
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by ?/o and the proportion of saturated stations by yx- Denoting p = /i(s — 
Sfce{o k} kyk)/\ these two quantities are given by: 

V°0>) = p K+i _ ! and V*(P) = P R -pkiZTi ( 7 ) 

One important problem is to study the effect of the number of bikes per 
station s on the performance of the system given by equation (7). Even if the 
fixed point of the equation (4) can not be written in close form as a function 
of s as soon as K exceeds 4, equation (7) provides an efficient way to achieve a 
performance study of the system by considering the parametric curve 

(- + X>Ao(p), yo(p) + y K (p)) (8) 

V ^ k=l ) { p >o} 

The use of this parametric curve allows us to study efficiently the performance 
of the system as a function of the number of bikes per station s. These results 
are mainly summarized in the next theorem and on Figure 1. 

Theorem 2. For the one cluster model, 

(i) the limiting proportion of problematic stations, y$ + 2/K> * s minimal when 
s = K/2 + X/fj, and the minimum is equal to 2/(K + 1). It goes to one 
when s goes to zero or infinity. 

(ii) As K grows, the performance around s = K/2 + X/fi becomes flatter and 
insensitive to X/p. 

Proof. Let <p(p) = y (p) + y K {p) = 1 - (p K - p)/(p K+1 - 1)) and s(p) = 
Ap//i + [I, Up]. Functions if and s are well defined on R+. The proportion of 
problematic stations as a function of s is given by ip = ipo s^ 1 . First prove that 
ip has a minimum at sq = s(l) which is 2j(K + 1). Differentiating function ip 
with respect to p gives 



p 2K -l + K(p K - 1 - p K+1 ) 

{p K + 1 - l) 2 



Differentiating the numerator and studying the variation, it holds that the nu- 
merator is strictly increasing on E+ and <p'(l) = 0. Thus tp is strictly decreas- 
ing on ]0, 1[, strictly increasing on ]1, +oo[ and has a minimum at 1 which is 
(f(l) = 2/(K+l). It shows the optimal number of bikes per station corresponds 
to p = 1, and thus s(l) = K/2 + X/p is the optimal proportion of bikes per 
station. This leads to a proportion of problematic stations of 2/(K + 1) and 
concludes the proof of (i). 

To prove (ii), we compute the second derivative of ip = (po s _1 at sq = s(l). 
A direct computation shows that ip"(so) — tp" (1) / s' (I) 2 . Replacing p by 1 + x. 
we derive the asymptotic expansion of (p at p = 1 

2 1 K(K - 1) . n9 

ip(p) = + — ^ Up - 1) 2 + o((p - i 3 ) 
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and therefore that (p"(l) = K(K—1)/3(K+1). Moreover, differentiating s gives 
that s'(l) = A//x + K 2 /12 + K/6 which leads to 

V (so) - 



3(K + l){\/n + K 2 /12 + K/6) 2 ' 
It allows to estimate the behavior around the minimum. First, 



K{K + 2f 



which means that ip can never be sharp for the range of values considered here, 
i.e. K > 20. The quantity ip"( s o) is increasing as a function of X/fi. If X/fi 
is negligible compared to K 2 /12, the behavior does not depend on X//I. More 
precisely, if K is large, for A//x fixed, ip"(so) ~ 48/K 3 . □ 

This theorem indicates us that even for a symmetric system for which the 
number of bikes per station is chosen knowing all parameters of the users, the 
proportion of problematic stations only decreases at rate 1/K. This is problem- 
atic for practical situations where because of space constraints and construction 
costs, a station cannot host more than 20 or 30 bikes. A system with 30 bikes 
per station would lead to a proportion of problematic stations of 2/31 rj 6.5%. 
Although it might be acceptable if this bike-sharing system is used once in a 
while, a probability of 6.5% of problematic stations is too large for a reliable 
daily mode of transportation. 




Figure 1: Proportion of problematic stations as a function of the number of 
bikes per station for various values of the size of stations K. On the cc-axis is 
the average number of bike per station s. For both scenarios, we plot \/fi = 1 
and X/fi = 10. 

The results of Theorem 2 are completed by Figure 1, that shows how the 
proportion of problematic stations evolves as a function of the number of bikes 
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per station for K = 30 and K = 100. For both scenarios, we plot the results for 
A//i = 1 and X/fi = 10. 

When the number of bikes that can be hosted per station is fixed to 30 
and A//i = 1, the performance is almost equally good from 10 to 20 bikes per 
station. For Velib, s = 20000/1500 ~ 13. However, as soon as the number of 
bikes per station is lower or bigger, the performance decreases a lot. When the 
size of the stations is K = 100, the performance is less sensitive to the number 
of bikes per station and is slightly better. As pointed out by Theorem 2, in that 
case the proportion of problematic stations is 2/(K + 1) 2%: Multiplying the 
maximum number of bikes than can be hosted at a station by 3 divides by 3 
the proportion of problematic stations. Having some stations designed to host 
up to 100 bikes is realistic for stations near subway for example but having all 
stations in a city with 100 spaces is very costly in terms of space and installation 
cost. In the following, all our simulations will be made for K = 30. 

When A//i = 10, the situation is similar to the one when X/fi = 1, with 
curves shifted to the right. The minimum number of problematic stations is 
the same, only the optimal fleet size is changed. This fact can be deduced from 
(8) where the term Xj \i only affects s(p) = X/pp + J2^P k Vo(p) an d not the 
proportion of problematic stations yo(p) + Uk{p)- 

These results suggest that without incentives for users to return the bike to 
a non-saturated station or any load-balancing mechanism, the implementation 
of a bike-sharing system will always observe a poor performance, even if the 
system is completely symmetric and there are no preferred areas. In a real 
system where some regions are more crowded than others {e.g. because of the 
trips from housing areas to working area), the situation can only be worse. 
In the following, we will examine mechanisms than improve dramatically the 
situation and that could be implemented at little cost. 

3.3 If people return the bike at a non saturated station 

Before studying incentives or regulation mechanism, we first study a variant of 
the model where users know which stations are empty or full. They arrive to a 
non-empty station and return their bike only to non-saturated stations. 

The dynamics of the system are slightly modified as follows. As before, there 
is a Poisson arrival process in the system with rate NX, but each arriving user 
chooses among the stations a non empty station at random. If there is no non 
empty station, he leaves the system. If the user manages to take a bike, then 
after a time exponentially distributed with parameter p, he arrives at a non 
saturated station chosen at random (i.e. with less than K bikes), returns his 
bike at this station and leaves the system. If there is no non saturated station, 
the user will try again after some random time, exponentially distributed with 
parameter p as in the previous model. 
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The transitions of (Y(t)) are now given by 



V -> V + ^( e k-i - efc) ^ XVk Nl k>0tyo<1 (9) 

iv l — 2/o 



y -> 2/ + -^( e fe+i - e fe) i /* 

iV 1 -2/K 



(10) 

V n=0 / 



These transitions are similar to the ones of the original system described in 
Section 2.2 and are transitions of the form y — > y + £/N which occurs at rate 
Nj3e(y) for a finite number of I g C. Let / be defined by f(y) = Yleec ^^(v) : 



^ f x f \ i , [Ay]), 



/(f) = Z^M i ( efc -i ~ e k)H>o H 1 ' (efe+i - efc)lfc <if 

As opposed to Section 2.2, the function f(y) — X^e/;^/^(y) is discontinuous 
when 2/0 goes to 1 and when yx goes to 1. Notice that Kurtz' results are limited 
to Lipschitz-continuous functions. Nevertheless for a extension to take into 
account discontinuous dynamics, see [5]. The price to pay for this generality in 
that the classical ODE y — f(y) has to be replaced by a differential inclusion 
y € F(y) where F(y) is the convex hull of ^vev(y){f{ z ) ■ z ^ V} where V(y) 
is the set of neighborhoods of y. In particular, for all the point for which / is 
continuous, which is the case when yo ^ 1 and yx ^ 1, F(y) is equal to f(y). 

The differential inclusion that corresponds to F has a unique solution which 
is an absolutely continuous function such that y = f(y) almost everywhere. 
When s < A//z, the situation is simple: In the limit, the solution of the dynam- 
ical system is a trajectory which is constant with yo(t) = 1 after some time. 
It means that, after some time, for a large system, almost all the stations are 
empty: a bike left at a station is immediately taken by someone else. The case 
s > K + X/fjb is similar, all the stations being saturated. If s > X/fi, it can be 
proved elementary, using ODE y = f(y) for k = that, for all t > 0, yo(t) < 1 
and using it for k = K that, if s < A/ fi + K then, for all t > 0, yic(t) < 1. This 
latter condition is assumed in the following (in practice s < K otherwise the 
total fleet size exceeds the total capacity of the stations). 

Assume in the following that A/ /i < s < K + X/fi. Equation y = f(y) 
can also be rewritten as equation (3) where L y is the infinitesimal generator of 
an M/M/l/K queue with arrival rate /x(s — [I,y])/(1 — 2/if) and service rate 
A/(l — yo). The probability measure y on {0, . . . , K} is an equilibrium point 
of (y(t)) if and only if, using equation (3), yL y = that is y is the invariant 
probability measure v p ( y ) of a M/M/l/K queue with arrival-to-service rate ratio 

p(y) = Ml - yo)(s - [I, y])/(A(l - y K )). 

As in the previous section, let p satisfies the fixed point equation 

A \-v p (K) n 
fi 1-Vp{0) 
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Using that p(l — v p (K))/(l — v p (<Sj) = 1, this equation is equivalent to 



[I 



• "pi 



a/a*. 



(ii) 



It is then straightforward that this equation has a unique solution. Therefore, 
the dynamical system has a unique equilibrium point which is the geometric 
distribution on {0, ... , K} with mean s — X/p, denoted by v B ^ x l^ . 

Theorem 3. Let assume that s > X/p and let h be the relative entropy defined 
in (6). Then: 

(i) y I— > h(y / v s ~ x > v-) is a Lyapunov function on {y £ P({0, . . . , K}), [J, y] = 
s - X/p}. 

(ii) The Dirac mass at the equilibrium point is the concentration point of the 
sequence of the invariant measures as the system becomes large. 



Proof. The proof is given in appendix A. 2. 



□ 



When studying the fixed point of the system, the main difference with the 
original model studied in §3.2 is that s is given by Equation (11). Therefore, 
we can plot easily the performance of the system by a parametric curve similar 
to (8) but replacing the first term in X/ up by X/p: 



X 
/' 



K 



kp k y (p) , 2/o (p) + Vk (p) 



(12) 



fc=i / {p>0} 

When s ^ [X/p, K + X/p], the proportion of problematic stations is one. 




without forcing 

forcing to non-empty/non saturated 




Figure 2: Proportion of problematic stations as a function of the number of bikes 
per station when we force people to go to a non-saturated stations compared to 
the proportion if we do not force people. Values for K = 30 and X/p = 1. 

Equation (12) shows that the proportion of problematic stations has a similar 
shape for this model and the one of the previous section. As before, the minimal 
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proportion is 2/(K+ 1) and is attained for s = K/2 + When s is not equal 
to K/2 + A//i, this proportion will be higher than in the classical model. This 
is illustrated on Figure 2 where the proportion of problematic stations for that 
model is compared with the one of the model studied in the previous section. 
This shows that although forcing people to go to a non-saturated or non-empty 
station reduces the unhappy users since everyone can take or leave a bike at 
anytime, it makes the system more congested and does not improve the overall 
performance. A user that is not aware of the state of the system will have more 
problems in this system than in the classical one. 

4 Incentives and the power of two choices 

In this section, we consider that when a user wants to return his bike somewhere, 
it indicates two stations and the bike-sharing system indicates him which is the 
least loaded of the two. The technique of choosing the least loaded among two 
servers has been proved to be a very useful load balancing strategy [10] and is 
often referred to as the power of two choices. We show that in the bike sharing 
system, it also helps improving performance dramatically, even if only a small 
percentage of users obey this rule. 

4.1 The two choices model 

We consider a symmetric model with N bikes and s bikes per station. As before, 
users arrive with rate A in each station and take a bike if the station is not empty. 
However, when a user chooses its destination, instead of picking one station, it 
picks two stations at random and goes to the most empty. 

Let, for k 6 {0, . . . , K}, Uk(t) be the proportion of stations with more than k 
bikes. In Section 2, we have seen that the state of the system can be described by 
a vector y = (jjk)ke{o...K} with y k = 1, where yk is the proportion of stations 
having k bikes. Therefore, the state of the system can also be described by the 
vector u = {uk)ke{Q,...K} with u non increasing and uq = 1 since there is a 
one-to-one mapping from the set of possible y to the set of possible u defined 
by, for k, < k < K. u k = J2i> k Vi an( ^ conversely, for k, < k < K, 

Vk = u k - u k+ il k<K . Note that J2ke{o,...,K} k Vk = X)f=i u k- 

There are two types of transitions for the Markov process. The first one 
is a user that takes a bike from a station having k bikes. This happens at 
rate NX(u k — u k +i). The second type of transition is when a user returning a 
bike at a station having k bikes. Since a user chooses the least loaded of two 
stations, this happens at rate — u k v )(sN — J2 k u k)- This shows that as 

in Section 2, as N grows large, the behavior of the system can be approximated 
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by the dynamics of the following ODE 



u o (t)=0; 



(13) 



K 



u k (t) = -A(u fc (t) - u k+1 (t)l k<K ) + M(«k-i(t) ~ - J2 Uk ^' ( 14 ) 



A more realistic model is to assume that only a fraction of users follow this 
rule while the other ones would go to one station of their choice. This would 
happen for example if a user obeying to the two-choices rule would get a reward: 
only a fraction of users would value this reward more than the cost of not going 
exactly where they want to go. Similar incentives exist in the Velib system in 
Paris (see footnote 2 on p. 20). To model this behavior, we assume that each 
user obeys to the 11 two choices rule" with probability r and otherwise chooses 
only one station and returns the bike at this one. 

The dynamics of the ODE can be written: For Uo(t) = 0; for k € {0, ... , K}, 

u k (t) = - X(u k (t) - u k+1 (t)l k<K ) 



Note that r = corresponds to the original ODE (2) while r = 1 corresponds 
to the case where everyone obeys the two choice rule, ODE (14). 

4.2 Equilibrium point of the dynamical system 

Let us denote p — p(s — Yl k =i u fc)M- An equilibrium point u of the ODE (15) 
satisfies the equations: For k € {1, ... , K}, 

u = 1; u k - u k+1 l k<K = p (r(u 2 k _ 1 - u 2 k ) + (1 - r)(it fe _i - u k )) . 

When p is fixed, we prove in Section A. 3 that this equation has a unique solution 
u{p). It is, as far as we know, an original result, and we write it in the context 
of queues with finite capacities. Moreover, we also show that u(p) is increasing 
in p. Let f(u) = p(s — J2 k =i U fe)/A. It should be clear that f(u) is decreasing 
in u. Therefore, this implies that p t— > f(u(p)) is decreasing. This implies that 
p = f(u(p)) has a unique solution on [0;oo). Therefore, the dynamical system 
(15) has a unique fixed point. 

4.3 Impact on the performance 

The proof of the uniqueness of the solution u(p), presented in Section A. 3, 
provides an efficient way to compute function of p. This shows that, 

if p is fixed, the number of bikes in the system is s(p) — Xp/p + Y] k _-, u k (p). 
This provides an easy way to plot performance indicators by using a parametric 



+ (r(t£_i(*) ~ "*(*)) + (1 - r)(« fc _i(t) - u fc (t))) 




(15) 
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curve of parameter p. These results are reported in Figure 3 and indicate that 
the performance of the system is radically improved, even if 20% of users obey 
the two choices rule. 




5 10 15 20 25 30 35 40 0.1 02 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Number of bikes per slation: s Number el bikes per station s 

(a) Proportion of problematic stations as a (b) Proportion of problematic stations as a 

function of the number of bikes per station function of the proportion of users obeying 

for r = 1 (everyone obey the rule) and r = 0. the two choices rule. 



Figure 3: Performance of the system in a two choice system: 3(a) when everyone 
obey the two choices rule and the number of bikes per station varies. On 3(b) 
when the number of bikes per station is 16 and r varies. 

On Figure 3(a) is reported the proportion of empty stations as well as the 
proportion of saturated stations as a function of the proportion of bikes per 
station when everyone follows the two-choices rule. Compared with the orig- 
inal case where no one follows the rule (see Figure 1), this situation has two 
advantages. Firstly the optimal performance of the system is much better than 
in the original system: Recall that K — 30. While in the original system, the 
proportion of problematic stations is at best around 7%, here the proportion of 
problematic stations can be as low as 10 -6 (this is the case when the number 
of bikes per station is approximately 23). The second advantage is that this 
curve is rather insensitive to variations in the number of bikes: the proportion 
of problematic stations is less than 10 -3 if s is between 10 and 27 bikes. 

On Figure 3(b), we fixed the number of bikes per station to 16 and let the 
proportion of users obeying the two choices rule vary from to 1. We choose 
16 since this is the optimal fleet size without incentives (see Section 3). This 
figure shows that the proportion of problematic stations diminishes rapidly as 
soon as the number of users obeying the rule grows. Moreover, the decrease of 
the log of this proportion is almost linear: if 25% more users obey the rule, the 
proportion of problematic stations is roughly divided by 10. 

These results suggest that incentive mechanisms of that type should be im- 
plemented in practice. Even if the advice is only followed by a few users, this 
leads to a dramatical improvement of the situation. This point will be further 
developed in Section 5.3 where we study the performance of incentives against 
the performance of regulation mechanisms. 
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5 Asymmetric model and regulations 



In practical situations, the system is not symmetric. This may happen if the city 
is composed of areas more occupied by housing while other areas concentrate 
more working places. In that case, depending on the time of the day, there are 
routes that are preferred to others: e.g. in the morning, people go from their 
home to their working place. This can also be the case if some stations have a 
higher elevation than others. In that case, many people want to take a bike to 
go down but very few are willing to take a bike and ride up. 

In this section, we consider a model with two different clusters, that can 
represent either housing and working area in the city or stations with various 
elevation. We first study the steady-state performance of the system in the 
absence of regulation, showing that the performance is very bad even if the 
system is slightly asymmetric. Then, we show that adding trucks that balance 
the number of bikes between stations improves the situation a lot, if the number 
of trucks is large enough. 

5.1 Model and performance in the absence of regulation 

In Section 3, we have seen that when the system is symmetric, the performance 
is optimal when the number of bikes per station is slightly more than half of the 
available spaces at the station. Here we show that in an asymmetric case, the 
number of bikes should be less than that but that, in any case, the performance 
is much worse than in the symmetric case. 

We consider a model with two clusters. People arrive in cluster 1 with rate 
Ai and in cluster 2 with rate A2. When a user takes a bike in area i, it chooses 
with probability one half to go to a station of area j 7^ i and one half to go to 
a station of area i. If a user arrives at a saturated station in area i, he goes 
to a station of the same area at random. The mean travel time between two 
stations is equal to = 1, the stations being in the same area or not. Again, 
the number of bikes K than can be hosted in a station is 30. 

As opposed to Section 3, there is no direct scheme for proving the unique- 
ness of the equilibrium point of the ODE representing the limiting dynamical 
system. However, this uniqueness is conjectured and confirmed by numerical 
experiments. The fixed point of the ODE can be easily computed by simulating 
the ODE. For each cluster, the proportions of stations to be either empty or 
saturated are reported in Figure 4 for two values of A1/A2, 1.3 and 2 (recall that 
M =l). 

In both cases, the proportion of saturated stations of the overloaded area 
starts growing when the number of bike per station is around 10 to 15 and goes 
to 1 very quickly after that threshold. On the other hand, the proportion of 
saturated stations of the underloaded area remains always small, except when 
the number of bike is extremely high: more than 40 bikes per stations with 30 
spots even for X1/X2 — 1.3. This is explained by the huge number of people 
that are in the overloaded area looking for an available station. 

We also computed the proportion of problematic stations for an asymmetric 
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Empty stations of cluster 1 Empty stations of cluster 1 

Saturated stations of cluster 1 " 09 Saturated stations of cluster 1 

Empty stations of cluster 2 Empty stations of cluster 2 

Saturated stations of cluster 2 " Saturated stations of cluster 2 




(a) aT/As =™1*3 (b) N l7/A7=2 



Figure 4: Performance of the system as a function of the number of bikes in a 
city with two areas that have different arrival rates. 

system with A1/A2 = 1.3 and 2. These results are not shown here due to space 
constraints and can be obtained directly from Figure 4 by computing the sum 
of the four curves. These has two implications on the design of an efficient bike- 
sharing system. The first one is that the optimal number of bikes per station 
is less than in the symmetric case: around 13 for A = 1.3 and around 9 for 
A1/A2 = 2. That suggests that in a real system where A1/A2 is unknown, the 
number of bikes should be a little less than half of the available spaces but 
not more. The other and probably main conclusion is that asymmetry without 
regulation implies a poor performance of the system. Even when the system is 
close to symmetry (A1/A2 = 1-3), if the number of bikes is optimal, a user going 
from an underloaded to an overloaded station has more than 25% chance of not 
finding a bike and more than 10% of finding a saturated station. 

5.2 Balancing the load of stations by trucks 

The problem of balancing the bikes in various areas in a city is one of the 
major issue for the people who conceives such bike-sharing systems. One simple 
solution is to use trucks that move bikes from saturated stations to empty ones. 
In this section, we study a simple load balancing strategy where trucks move 
at random. We show that when the rate at which the trucks visit station is of 
order 10% of the arrival rate of users, the performance is improved dramatically. 

We consider the same model as in the previous paragraph adding trucks 
that move at random. We consider a fleet of 7 trucks per station. At rate 1, 
each truck chooses two stations at random and equalizes the number of bikes 
in the two stations. We assume that the average time for a truck between two 
equalizations is the same as the average time of a user to go from one station 
to another one which is /i^ 1 = 1. This is due to the fact that loading and 
unloading bikes take time. This means that two stations are equalized at rate 
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N-f 1 . For simplicity, we assume that the bikes are moved instantaneously from 
one station to another. This means that, at rate N-f, two stations in states (c, i) 
and (c',j) become in states (c, [(i + j)/2j) and (c', \(i + j)/2}). 




Figure 5: Performance of the system as a function of the ratio of arrival rates 
in the two clusters A1/A2 for various values of the number of trucks per station 
7- 

On Figure 5 are reported the sum of the proportion of problematic stations 
for the two clusters for various values of the ratio of arrival rates in the two 
clusters A1/A2 and for a number of trucks per station of 7 = 0.01 and 0.1. As 
mentioned in the previous section, without regulation, the performance of the 
system decreases dramatically as soon as A1/A2 grows. Figure 5 shows that 
for moderate values of A1/A2, 7 = 0.01 trucks per station does improve the 
performance. However, this improvement saturates quickly for higher values 
of A1/A2. When the number of trucks per station is higher - 7 = 0.1 - the 
performance of the system degrades very slowly with A1/A2, showing that 7 = 
0.1 really improves the performance. 

5.3 Comparison of incentives and regulation 

In Section 4 dealing with a symmetric system, the situation is greatly improved 
if the users return their bike at the less loaded of two stations chooesen at 
random. Here, we study how this strategy performs compared to the regulation 
mechanism presented in this part, both for a symmetric model (Figure 6(a)) 
and an asymmetry model (Figure 6(b)). 

On Figure 6(a), we compare the situation when all the users obey the two 
choice rule with a situation where 7 = 0.01 or 7 = 0.1 trucks per station 
balance the bikes between stations. Whereas the performance for 7 = .1 is 
quite good, it is still far from the performance of the two choice rule. However, 

1 For this model, we assume that there are 7 trucks per station moving at rate 1. If we 
consider that the moving rate of trucks is 8 and the proportion of trucks per station is 7/(5, 
this leads to the exact same behavior: two stations are equalized at rate JV7. 
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Figure 6: Performance of the incentive mechanism of Section 4 and the regu- 
lation mechanism for the symmetric case (one cluster) and an asymmetry case 
with two clusters. 



in an asymmetry model, the results are quite different and the two choice rule is 
not sufficient to insure a good performance. Figure 6(b) shows the proportion 
of problematic stations for the asymmetric model with two clusters presented 
in §5.1 with A1/A2 = 2 for various strategies. We compare four situations: the 
original model without regulation or incentive; two numbers of trucks per station 
(7 = .01 and 7 = .1); a situation where all the users obey the two choices rule; 
a situation where we combine regulation with 7 = .01 and where all the users 
obey the two choices rule. By two choice rule, we mean that a user chooses to 
return its bike to the least loaded station among two in the same cluster. 

The case where all the users obey the two choices rule improves the situation 
with or without regulation. However, for a number of trucks of 7 = .01 or lower, 
the performance is not good even in that case. The only configuration that 
provides a good performance is when the regulation rate is high. 

These results show a mechanism to balance the number of bikes between 
clusters is necessary to achieve a good performance. This mechanism can be a 
regulation mechanism as we studied here but it can also be incentives: users 
get a reward if they return their bikes in a station with a high elevation 2 . But 
the later, very interesting for a symmetric model, performs badly in presence of 
asymmetry. 



6 Conclusion and future work 

The basic large networks modelling bike sharing systems and taking into account 
the finite station capacities have been studied via a simple symmetric model for 
a mean field point of view, providing the limiting steady-state measure. The 
limit has a probabilistic interpretation which provides analytical expressions for 

2 Such a mechanism is implemented in Velib in Paris with the "V+" concept. V+ gives to 
the user 15 min to be used for a next ride if a bike is left at one of 100 uphill stations. 
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the performance of the system. The conclusion is that an optimal fleet size is 
highlight and performance is quite robust around this value for the range of 
values of station capacities considered. But the performance is not excellent. 

For that purpose, an incentive to return bikes to the least loaded station 
among two improves dramatically the performance, even if a small proportion 
of users accept to do this. Regulation could help if the number of trucks is 
sufficient. An asymmetric example is proposed via clusters, proving the dire 
need to improve the poor performance of the system in this case. Although 
clusters make the limiting dynamical system analytically quite complicated, 
they allow a rather simple numerical evaluation of the system. 

The partitioning of stations into clusters provide insights on the effect of 
asymmetries in the system, such as overloaded and underloaded areas. However, 
a better understanding of the geometry is necessary. A user finding a full station 
is more likely to perform a random local search than going at another station at 
random. A similar situation arises for the incentive to choice the least loaded 
station among two. A more realistic model would be to consider a choice between 
two neighboring stations. Building a tractable model that includes geometry 
and heterogeneity is for the moment an open problem. This is an important 
question that should be addressed in future work. 
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A Appendix 



A.l Proof of Theorem 1 

The proof mimics the proof by of Proposition 4.1 of [12]. The main difference 
is that the invariant measure v p associated to L is geometric but not Poisson, 
but it does not change the arguments. Using that v p (n) — p n /Z(p), one has: 



h {y/ v P (y))= 2/™ lo s 



y n Z{p(y)) 



p(y) n 

n£{0,...,K} IKy> 

Vn log(y„) + log Z{p{y)) - [I, y] \og{p(y)). 

ne{0,...,K} 



Plugging in g gives 



g(y)= Y yn l °&(yn) - [i,y]iog(p(y)) + i>([i,y]). 

ne{0,...,K} 

Thus, using that, for each n, ^—[1, y] — n, that p{y) — (p([I, y\) and that tp is 
a primitive of x — > Xip'(x)/ip(x), this leads to: 

^-9(V) = log (^) + 1 = log + log Z(p(y)) + 1 

"Vn \PKVrJ \ v p{y){ n )J 

Let q y {m,n) = v p ^(m)L y (m,n). By reversibility of generator L y with respect 
to distribution v p ty), <ly(m,n) is non-negative for m ^ n and symmetric in 
(m, n). Using the symmetry of q y (.,.) and the fact that L y (m,n) = 0, a 
direct computation shows that, for all u £ Kl°^- ^1, 

V q y (m,n)( — ^— -^r ) (u m - u n ) 

= -2 Y y m Ly(m,n)u n 

n,me{0,...,K} 
= -2{yL y ,u) . 

As a consequence, 
{yL y ,Vg(y)) 

1 V - * / \ ( Vm Vn \ f, Vva , Vn 

= — — > QyiTfl, n) ~, r -j-^ log -. — log ^— ■ 

2 {m,n)i{£,...,Ky \"p(y)( m ) »P(y)( n )J \ v P (y)K m ) " P (y)W 
= - o 2^ <ly( m i n ) 7 — r — 1 log ~T^- 

2 (m,n) S R...,/f} 2 V P{vm \ V P{y)\ m ) Vn J V^p{y) M 

Each term of the sum is either equal to when m — n or of the form C(x— 1) log x 
which is non-negative for each x > as C > 0. This shows that: 

J t 9(v(t)) = (y(t),Vg(y(t))) = (y(t)L y{t) ,Vg(y(t))) < 0. 
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Moreover, the derivative of g is equal to when for all n,m: y m u p / y \(ri) = 
Un v p(y){ r n)- Thus there exists a constant C such that y n = Cv p ^ y ){n). Since 
12 Vn = J2 v p(y)( n ) = 1> Vn = v p (y){n) for all n. □ 

A. 2 Proof of Theorem 3 

The proof of (i) is follows the same lines as Tibi [12, Proposition 4.3]. It will 
be omitted. 

The proof of (ii) can be deduced from (i), mimicking the proof in Tibi [12, 
Proposition 4.4]. Function h(y/v p t y -j) is not a Lyapunov function on V({0, . . . , K}) 
but only on V({0, . . . , K}) s _ x/p = {y e V({0, ... , K}), [I, y] = a- A/a*}, by (i). 
Nevertheless the quantity l(y) = ([J, y] — s + A//i) 2 decreases along the flow. To 
see this, 

K-l 

< yL v> Vf(y) > = - ^2 (yn+iL y {n + l,n) - y n L y (n,n + l)) 

71=0 

= - K f(^±l-y n ^ {s - {I , y]) ) 

= -(A -[/,«])) 
= -p{[I,v] -X/fJt + s) 

where f(y) = [I,y] = J2ne{o,...,K} n Vn- Thus, < yL v ,Vl(y) >= 2([I,y] -s + 
r) < yL y ,Vf(y) >= -2(j,([I,y] - s + \/n) 2 < 0. We conclude as in [12], that 
any weak limit of the sequence of invariant measure is supported by the set 
V{{Q,...,K}) S _ X/}1 = {y e V({0,...,K}),l(y) = 0}. As y ^ h{y/u p(y) is 
a Lyapunov function for the system on V({0, . . . , K}) s _ x / pi it is proved that 
the Dirac mass at v s ~ x /l l is the only weak limit for the sequence of invariant 
measures. It proves convergence by compactness of ^({0, . . . , K}). □ 

A. 3 Equilibrium point in the supermarket model with fi- 
nite capacity 

The model is a queueing system of N queues with finite capacity K where the 
arrival process to the system is Poisson with rate AN and the arrival customer 
either chooses with probability r d queues at random and joins the shortest 
queue, ties being solved at random, or joins with probability 1 — r a queue at 
random. When r — 1, the model is the classical model where the customer 
joins the shortest queue. But it is not much difficult to deal with the general 
case where r € [0,1]. The following result claims the uniqueness of a fixed 
equilibrium point for the limiting dynamical system when N gets large. 

Theorem 4. For each p > 0, the system 

u = 1, u k - u k +il{k<i<} = P (rfak-i _ u t) + (1 - r)(uk-i - u k )) (1 < k < K) 

(16) 
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has a unique solution in U = {u = (uk)o<k<K € + , non increasing 7 uo — 
1}. 

Proof. Let uk+i be defined as follows: 

Uk+i = u K + p {r(u d K - uj i _ 1 ) + (1 - r)(u K - u K -\)) ■ (17) 

The key point is to notice that equations (16) and (17) are equivalent to 

u = 1, u k+1 = Ul + p(r(u d k - 1) + (1 - r)(u fc - 1)) (0 < k < K). (18) 

It is proved straightforwardly by induction on k. Then, let u\ — x € [0, 1] 
and define, for each k, u k as a function of x. It is thus sufficient to prove that 
uk+i(x) = has a root in [0, 1]. Indeed, such a root define a u £ U, solution 
of (16). Moreover the solution of (16) is unique. If for x and y with x < y, 
u(x) and u(y) are solutions of (16), then, straightforwardly with equation (18), 
Uk(x) < u k {y) for k < K + 1. It contradicts uk+i{x) — UK+i{y) = 0. 

For that, assume xi — and let us prove by induction that, for k £ 
{1, . . . ,K + 1}, x i — ^ u k {x) is strictly increasing on [2^-1, 1], < at x k -i, such 
that Mfe(l) = 1 > 0, that there exists x k unique root of u k (x) = on 1], 
and thus u k > on Jx^, 1] and x k < p. 

For k = 1, the assertion is true. Suppose the assertion true for k < K. Thus, 
for x € [x k , 1], by equation (18), 

u fc+ i(aO = x + p (r(u k (x) d - 1) + (1 - r)(u k (x) - 1)) . (19) 

and 

Ufc+iO") = 1 + P (rd« fc (a?) , '- 1 «;(a:) + (1 - r)<(x)) . 

with, by induction assumption, u k > on ]x k , 1] and u' k > on ]x k , 1] because 
u k is strictly increasing on ]xk-i, 1]- Thus u^+i is strictly increasing on [x k , 1]. 
Moreover u/-_|_i(l) = 1 by equation (19) and induction assumption. By equation 
(19) and by definition of x kl u k +\{x k ) — x k — p < because x k < p. Thus there 
is a unique x k+ \ G]a;fc,l[ such that u k +i{x k +\) — 0. Moreover, u k+ i > on 
[x fc+ i,l]. By equation (19), X k+1 = p(r(l-uf(x k+1 )) + (l-r)(l-u k (x k+1 ))) < 
p because u k (x k+ i) > due to x k+ \ €]x k _i, 1[. Thus xk+i is a root of % + i 
in [0,1]. It concludes the proof. □ 
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